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Renewable energies are in the news these days and are seen as a solution to 
endorse energy independence and diminish greenhouse releases. One 
possible new renewable source that is emerging as a promising technology is 
the fuel cell. The ‘F.C’ is an electrochemical system that transfers the 
chemical power of a redox response into electrical energy with simultaneous 
production of water and heat. In fact, there are numerous kinds of fuel cells; 
in our work we are interested in studying the solid oxide fuel cellSOFC. In 
the context of this work, a modeling tool has been implemented SOFCs to 
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cen the analogy between the electrical, thermal and chemical domains; this way 
Electrochemical of proceeding constitutes a simple, evolutionary and efficient tool. Using 
Fuel cell this modeling, a simulation was carried out in order to obtain the dissimilar 
Modeling characteristics and the impact of the studied parameters on the performance 
SOFC of the SOFC fuel cell. 
Solid oxide P . . 
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1. INTRODUCTION 

In recent decades, solid oxide fuel cell SOFCs cells have attracted a lot of attention, owing to their 
potential uses as stationary power producers and mobility (ground, marine, air). The great energy conversion 
efficiency and minimal hazardous emission levels of SOFCs make them appealing (only the CO2 released by 
the hydrogen production process is a concern). Modularity, fuel flexibility, and minimal noise are among the 
other benefits [1]-[3]. Furthermore, high operating temperatures afford good features, such as potential utilize 
of SOFCs in extremely proficient cogeneration applications. Solid oxide fuel cells are also appropriate for 
internally reforming the fuel (e.g.natural-gas, propane, methanol, gasoline, and diesel), as a result making it 
doable to stay away from the acceptance of both highly sophisticated, expensive outside reformer and to 
remove the complexities from fuel storage [4]. The giant confronts to promote the dispersion of SOFC-based 
energy conversion systems are mostly connected to manufacture costs and permanence. The 
accomplishments of these objectives shall definitely contribute to upholding the technology and lastly 
beginning a mass production stage. Above and beyond expenses and performance, long-standing steadiness 
is a crucial condition for the marketable application of the SOFC technology [5], [6]. 

The hydrogen vector, namely fuel cells and electrolysis (to make hydrogen), also offers promising 
potential for these large-scale storage applications, where conventional batteries play a key role. As a result, 
the storage solution will be multi-system, and the electrical study of the association of many sources of 
energy - hybridization - will take on new significance in this context. Solid oxide fuel cells (SOFC) are 
electrochemical devices that allow chemical energy to be converted into electrical energy. The premise is 
straightforward: water and electricity are produced from a fuel like hydrogen and an oxidant like oxygen. 
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Apart from the fact that these systems are environmentally friendly since they do not emit greenhouse gases, 
they also show the potential of hydrogen as an energy vector for both the production of electricity and the 
storage of energy in chemical form. The applications primarily pertain to decentralized energy generation for 
residential and urban regions with electrical powers ranging from a few kW to a few hundred kW over 
hundreds of hours [7], [8] 


2. OPERATING MODEL OF SOFC CELL 

Cell is a set o of an anode, an electrolyte, a cathode and the interconnections. The anode is supplied 
with hydrogen or gas mixture (fuel) and the cathode is supplied with air (oxidizer). The operating principle of 
the SOFC is hinged on the following mechanism: oxygen is dissociated at the cathode into O2-, then 
the cation migrates through the ionic conductive electrolyte at high temperature (from 700 °C) and will 
combine at the anode with hydrogen to form water and release electrons that go via an outsider electrical 
circuit [8]-[10]. The main characteristic of SOFC is therefore their operating temperature (600 to 1,000 °C), 
which is mandatory to obtain a adequate ionic dynamism of the ceramic electrolyte. This temperature has a 
double advantage. Firstly, it allows the direct use of hydrocarbons, which can be easily reformed without the 
need for a noble metal catalyst. Secondly, it produces a high level of heat that can be easily used in 
cogeneration, with an efficiency of up to 80%. 

But it also has a disadvantage; the heating is long and complicates any use with short and repetitive 
cycles (as in the case of transport). For these reasons, the technology lends it itself particularly well to 
decentralized electricity production and cogeneration (areas covering powers ranging from 1 KW to several 
tens of MW) [11]-[13]. Thanks to its high efficiency and its potential capacity to operate directly with liquid 
hydrocarbons, it will also find an outlet in naval propulsion. The development of this type of cell, because of 
its high operating warmth and the resolution of thermo-mechanical problems of material resistance, is quite 
complex. One of the particularities of the SOFC is its solid electrolyte, usually zirconium (Zr2) doped with 
a mole of 8 to 10% ytterbium (Y3+),which acts as a conductor for the oxygen anion (O2-).The following 
Figure | shows the operating principle of the SOFC fuel cell. 
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Figure 1. Picture of SOFCs 


3. POLARIZATIONS 

Polarizations, or overvoltage’s, are voltage losses as consequence of drawbacks in thestructure of 
the F.C polyanionic cellulose (PAC) and materials microstructure. The ohmic confrontation of ‘oxyg’ ions 
passing through the electrolyte (i.R Q), electro_chemicals activated hurdles at the ‘ano’ and ‘cath’, and 
terminaly focus polarizations in view the incapacityty of gas’s to distribute at high velocities through the 
porous anode and cathode [10]. The cell voltage can be calculated employing in (1): 


Volt = Ey — R * i — Ncatnode — Nanode a) 


In SOFC, it is frequently essential to spotlight on Ohmic and attentiveness biases because high 
operating heats background little’s activation polarizations. Nonetheless, like the lowest boundary of the 
S.O.F.C functioning hotness is came close to (~873,15 kelvin), these polarizations turn out to be significant. 
Research by Huo [14] The above equation is utilized to find out the S.O.F.C voltage. This method showed 
results which gives aexcellent accordwith scrupulous investigational Information and deprived agreement for 
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experimental working parameters other than the originals. Additionally, most of the equations employed need 
the addition of multiple factors that are hard or unfeasible in order to identify. This that any process of 
optimizing the S.O.F.C operating facto rand the choosingvery complicated geometry configuration. Because 
these circumstances, a few other equations have been suggested [15]. 


Emax—tmax-Nf 11 
Esorce = -rin n (2) 
r(1-nf)+1 
Where: 
E_sorc : Cell voltage 
E_ max : The Nernst equation predicts the maximum voltage 


İ max : Current density at its maximum (for a given fuel flow) 
nf : Factor of fuel usage [15], [16] 

ra : The electrolyte's particular ionic resistance 

r2 : The electrolyte's particular electrical resistance 


4. THE LOSSES SOFC MODEL 

There are three main forms of polarization overvoltage’s: activation, ohmic and concentration. A 
constant offset also contributes to the total polarization, resulting from small overvoltages like contact 
resistor, insidecurrenty and leakage. According to [17], the shift was approximated to be 0.07 volts. the 
additions of the different polarizations result in a voltage sag from the optimal Nernst potential to the actual 
working value. As the interconnector and the electrodes are iso-polar, the array voltage is constant over the 
cell and can be expressed in (3): 


Vsorc = Ebers z Viet = Vonm or Véon a Vå ffset (3) 
The total energy derived from the S.O.F.C is determined as shown in (4): 


Ecet = Vsorc-!sorc (4) 


4.1. Activation polarization 
‘Activ? polarization means the power obstaclebeing diminished in order to  enablethe 
electro_chemical feedback taking partat the interface of the electrodes [18]. This quantity of power 
unavoidablyleads to a noteworthyovervoltage, that is typicallymodelled by the non-linear equation 
recognized as butler-wolmer formula [19]: 
RTE 


Vice = 2 (sinh (L) (5) 


a(T!).F 2Jo(T!) 
where æ is the charge transmit coefficient and J is the switch present mass. 


4.2. ‘Ohm’ polarization 

‘Ohm’ polarization primarily relies on the ionic momentum of electrodes and the electronic 
momentum of an electrolyte. Such over voltages are assessed finalizing the content of every part of the 
S.O.F.C, as (6): 


Vonme = Fan T (6) 


Vonm = È Vénmx k = [ano, cath, el] 
The conductances are approximated with, suggestions in the review of literature [17] for 2nd edition ceramic 


'S.O.F.C' stacks, in which the ‘ano’, ‘cath’ and ‘el’ cover materials are nickel cermets, strontium-doped 
lanthanum manganite and yttrium oxide cured zirconia, correspondingly: 


Oan = 10? 
Oca (T;) = C4 (Ti)? = C2 (TÌ)? = C3 (7) 


Go(T;) =c (Ti — 273.15)? + c_s(Ti — 273.15) + c6 
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4.3. Concentration polarization 

As the relative stresses of hydrogen and oxygen at the ‘ano’ and ‘cath’, respectively, drop as fuel is 
depleted. The rate of depletion is determined by the cell's average current density. When a result, as a current 
thickness rises, the relative stresses fall, and a deficiency number re-agents is eventually carried to the 
electrodes. Until the voltage is decreased to [20], [21], considerable losses occur. ‘Ano’ and ‘cath’ restrict 
currents are thevalues at which this phenomena occurs. Concentration polarization is the term for the voltage 
loss that occurs at high current densities and can be calculated as (8). 


Vim =~ 22. [tin (1-4) +n (1-4) - in(1 + Fig )| (8) 


2.F Jas PHoo-Scs 


The ‘ano’ and ‘cath’ limiting currents are determined as a function of the species scattering coefficients, 
using the formula given in [22]. 


5. MODELING OF THE SOFC 

The purpose of this paper is to give the model adopted to account for the physical phenomena in the 
cell. The first part of the paper is dedicated to the modeling principles applied to the chemical, thermal and 
electrical domains. The second part concerns the application of the modeling to the S.O.F.C Cell. The 
chemical § thermal modeling of the ‘F.C’ is based on their formal analogy with the electrical domain. The 
temperature is similar to the voltage and the heat flow is similar to the current. Two heat sources within the 
cell are considered: endothermic or exothermic heat sources produced by chemical reactions and Ohmic heat 
sources. Furthermore, pressure is similar to voltage and flow is similar to current. The application of this 
analogy offers the possibility of representing physical laws within a system with electrical circuits. The use 
of such a representation allows the use of different techniques in the electrical domain: putting chemical and 
thermal models in series or in parallel, realization of the regulation of chemical or thermal quantities with the 
help of methods used in automation (ex: pole compensation), ease and flexibility of integration of developed 
models in a global model (battery coupling or battery associated with the converter) The electrical model 
translates the electrochemical phenomena of each cell. The electrical model of the SOFC is based on the 
reversible Nernst voltage and the different voltage drops. 


5.1. Electro-chemical system 

We will choose a static SOFC system based on thermodynamics and electrochemical reactions to 
determine the electrical characteristics of the fuel cell.In a second step we vary the parameters and observe 
the variation [23]: 


(AG = AH—TAs) (9) 
(AH = AU+ APV) (10) 
(AU =w+Q) (11) 


Av: Growth in size, AU: Internal energy changes, œ: labour, ọ: thermal. 


(AG = w + ọ + PAV — TAs) (12) 
(Q = TAs) (13) 
(AG = w+ PAV) (14) 
(-PAV), thus: 

(AG = —nl'e.—PAV + PAV) (15) 


As a result, it was simply: 
(AG = —-nfe.) (16) 


With: 
e: electromotive strength of the stack, T: Constant of Faraday, n: number of electrons transported. 
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AG, = Gmo — Gm, — 5 Go, (17) 
1 
AG, = G, — NRT log ae (18) 
Po,» Pu,» Poo 
i 
E = £9 4 Zog te (19) 


nF PH 20 


5.2. Real potential of the stack 

The survoltages happened because of the shedding of non-reversible load, also named, which are 
polarization of ‘activ’, ‘Ohm’ polarization § ‘conc’. High-temperature fuel cells are able to generate enough 
heat to power the vaporiser, to power the steam generator. This will produce a hydrogen-rich refonnate which 
is then which is then supplied to the cell. 


V= E- losses (20) 


5.3. Capacity and quality requirements 
When the temperature of the battery drops, the variance of free energy increases, hence the standard 
voltage of a stack is in relation to the change in its standard free energy [24]. 
For water produced in a liquid state at 25 °C, it is: 
AGr = -237,190 joul/Mol. 
AHr = -285,840 joul/Mol. 
The Rated stack Voltage is: 


Ey = -22 = 122.9. V 

AGr = -228,590 j/Mol. 

AHr = -241,830/mol => E0 = 11.8 V. 
However at 1000 °C [25]. 


1000 
AHiooosc = Ahzsec + fos ACpdt (21) 


1 


000 C 
AS; o00°¢ = S25% + J, dt (22) 


5 
AGyo00°¢ = AH1000°c — TAS 000 °c (23) 
In this scenario, the start signal voltage is: Vo = 0.935 volt. 


5.4. Performance of the cell 

The chemical reaction is followed by a change in entropy, which leaves a free enthalpy, which may 
be transformed into work, in this case electrical energy, with the remainder being converted into heat during 
the chemical reaction in the stack. The report of the free enthalpy on the standard enthalpy of the reaction of 
water formation corresponds to the performance thermodynamics for a reversible transformation: 


AG° 

AH max me 
This performance ranges from 83 percent for liquid water creation to 95 percent for water formation in 
gaseous forms under typical conditions, but above 10° °C, the theoretical performance is: 


Ltn = 0.7337% 


a. The voltaic performance [24]-[26]: 
Of the in (24) one can define a voltaic performance to performance by voltage: 


vD 
Ue ==> 


(25) 


£0 
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b. The electrical performance: 
It is motivating to describe the notion of electrical performance: 


AG °Ve(1) 


[etectrique =TE AH!E® max (26) 
nF max 
Ve) 
[eectrique = or (27) 


It is then possible to directly measure this electrical performance experimentally using a voltage 
measurement, or to plot the electrical performance to T Constant by drawing the characteristic V(I) returns. 
Because this performance is highly dependent on temperature, it is necessary to recalculate Emax for each 
temperature [27]: 


Ve) Ve) 
Tronctionnemnt = aie = Ties (28) 


6. SIMULATION OF THE SOFC 

Simulation and mathematical models shall undoubtedly help in the development of a variety of 
power generators. But they are even more important in the development of fuel cells. Simulation of these 
cells allows you to focus on experimental research and improve the accuracy of interpolation and 
extrapolation of the results. In addition, mathematical models serve as a valuable tool for the design and 
optimization of fuel cell systems. Simulations are run on different systems with different software packages 
[28]-[30]. The following Table | shows the simulation parameters. 


Table 1.The SOFC model values [30] 


Settings Data Available units 

Nominal power 15,000 Watt 
Number of Stack 384.00 / 
‘Ano’ forces 2,280.006 Torr 
‘Cath’ forces 2,280.006 Torr 
fuel airflow 0.0012 Mol/s 
combustible Water rate 0.00001 Mol/s 
combustible Hydrogen rate 0.00009 Mol/s 
T fuel 799.85 °C 
T air 799.85 °C 
Current Density 60.00 A/cm? 
Vout 440.00 Volt 


7. RESULTS AND DISCUSSION 

In our study by using a simulation model with block diagrams under simulink Matlab to obtain the 
electrical and electrochemical characteristics of the SOFC type fuel cell, and then we varied the parameters 
of the effect of temperature and the effect of gas pressure, the results are the following: 

In order to verify the performance of the fuel cell generation system under normal operating 
conditions, the previously established model was simulated in the MATLAB/Simulink software environment. 
MATLAB/Simulink software environment. For this simulation, the model was planned With the parameter 
values extracted from the literature together with experimental data, the results are presented in the following 
and the interpretation is just after the presentation of the figures. 


7.1. Polarization curve 

A polarization curve for the SOFC fuel cell is shown in Figures 2 and 3. The characteristic begins 
with a no-load voltage at zero current, as can be shown. A dramatic reduction in voltage is noted when the 
cell current increases. The loss of activation voltage and parasites produced by the pace of the 
electrochemical process create this reduction in cell voltage. The non-linear activation losses contribute to the 
V(i) curve's downward tendency at low current densities. The activation losses remain the most significant in 
absolute terms. The Ohmic loss generated by the electrical resistances of the cell causes an almost linear drop 
in cell voltage and an increase in current above specific values [31]-[33]. 
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7.2. Temperature effect 

In this part we have varied the temperature for several values (T=600 °C, 700 °C, 800 °C, 900 °C, 
1000 °C) and we fix the partial pressure. Figures 4 and 5 show figures at various temperatures typical of 
SOFCs. The activation voltage drop dominates the voltage drop in the low current region (Figure 6). As the 
load current increases, the resistance voltage drop increases rapidly and is the main cause of the SOFC 
voltage drop (Figure 7). At certain load currents, the centralized voltage drop in the SOFC reduces the output 
voltage of the fuel cell Figure 8. Figure 4 also demonstrates the effect of temperature on SOFC-VI 
characteristics. The SOFC output voltage is high at low temperatures in the low current region and high at 
high temperatures in the high current region. From this figure, we can see that rising the cell temperature has 
a beneficial effect on the cell voltage. This is explained by the high ionic conductivity of the cells 
(electrolytes and electrodes) at high temperatures (such as T=900 °C). T=1000 °C) Figure 4. Therefore, it is 
needed to raise the temperature to ensure proper cell function [34], [35]. Tables 2 show the power values 
according to the influence of temperature: 
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Figure 8. Concentration losses characteristics for differenttemperaturesvalues 


Table 2. Influence of temperature on power 
TCC) Power (Watt) Ohmic Losses (Volt) Activation Losses (Volt) Concentration Losses (Volt) 


600 15,120 10.38 0.6345 0.001938 
700 15,25 12.58 / / 
800 15,32 15.8 / / 
900 15,35 20.8 / / 
1000 15,22 29.15 / / 


7.3. Pressure effect 
a. Variation of H2 pressure 

In this case the temperature is fixed at T=1000 °C and the partial pressure of the H2 gas is varied for 
the following values (PH2=0.01 Atm, 0.02A tm, 0.03 Atm) the results are the following. The polarization, 
activation, Ohmic resistance, and concentration-to-voltage curves at various pressures at 800 °C are shown in 
Figures 9, 10, 11, 12, and 13, respectively. In the event of a flooding failure, the pressure effect corresponds 
to the hydrogen and oxygen measured at the fuel cell inlet. 
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Figure 11. Activation losses characteristics for Figure 12. Ohmic losses characteristics for different 
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Figure 13. Concentration losses characteristicfor different pressures H2 values 


b. Variation of O2 pressure 

In this case the temperature is fixed at T=1000 °C and the partial pressure of the O2 gas is varied for 
the following values (PO2=0.01 Atm, 0.02 Atm, 0.03 Atm) the results are the following. As shown in Figures 
14, 15, 16, 17 and 18, as the pressure of oxygen increases, the activation loss decreases and the increasing 


current density decreases. In addition, increasing pressure reduces concentration loss and increases current 
density. 


=-<--- PO2,=0.01 Atm 
«E PO2,=0.02 Atm [T 14000 }- 
--#-- PO2,=0.03 Atm || 


--4-- PO2,=0.01 Atm 
+ PO2,=0.02 Atm 
=~ PO2,=0.03 Atm 


Voltage in Volt 
Power in Watt 
g 
s 


L i i $ | 
0 10 20 30 40 50 60 
Current Density (A/cm2) 


I | I i i 
0 10 20 30 40 50 60 
Current Density (A/cm2) 


Figure 14. Current-voltage characteristics for different Figure 15. Current-power characteristics for 
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Figure 16. Activation losses characteristics for different Figure 17. Ohmic losses characteristics for 
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c. Variation of H20 pressure 

In this case the temperature is fixed at T=1000°C and the partial pressure of the H2O gas is varied for 
the following values (PH20=0.1 Atm, 0.2 Atm, 0.3 Atm) the results are the following. As shown in Figures 19, 
20, 21, 22, and 23. In general, increasing the inlet pressure increases the fuel cell voltage and power density, 
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as shown in Figure 19. In normal operation of the cell, this pressure is maintained at 1 atmosphere. During 
the flood condition, the pressure at the inlet will automatically increase, causing cell failure. Therefore, 
reliable operation of the cell requires efficient monitoring of the cell. Tables 2 to 5 show the power values 


according to the influence of the pressure for H2O, O2, and H2. From the table, we can see that the value of 


power rises with the increasing temperature, and the parameter affecting this increase is Ohmic loss. After 
changing the pressure, you can see that the ohm loss and activation loss, and the concentration remain 
constant. The result was interpreted after the illustration of all figures. 
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Figure 18. Concentration losses characteristics for 
different pressures O2 values 
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Figure 20. Current-power characteristics for different 
pressures H20 values 
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Table 3. Influence of pressures (H2) on power 
Pm (Atm) Power (Watt) | Ohmic Losses (Volt) Activation Losses (Volt) Concentration Losses (Volt) 


0.01 21,150 10.38 0.6345 0.004275 
0.02 21,410 / / / 
0.03 21,630 / / / 


Table 4. Influence of pressures (02) on power 
Poo (Atm) Power (Watt) | Ohmic Losses (Volt) Activation Losses (Volt) Concentration Losses (Volt) 


0.01 17,540 10.38 0.6344 0.004275 
0.02 17,340 / / / 
0.03 17,120 / / / 


Table 5. Influence of pressures (H20) on power 
Pmo (Atm) Power (Watt) | Ohmic Losses (Volt) Activation Losses (Volt) Concentration Losses (Volt) 


0.1 14,790 10.38 0.6344 0 .004275 
0.2 14,600 / / / 
0.3 14,460 / / / 


8. CONCLUSION 

In the context of this work, a macroscopic modelling tool has been implemented thanks to the 
analogy between the electrical and thermal domains. This method constitutes a simple, scalable and efficient 
modelling tool. Once the component library has been completed, we will have a complete tool for studying 
and optimising the architecture of fuel cell systems, taking into account the different energy flows: electrical 
and thermal. Finally, we look at cogeneration systems using solid oxide fuel cells (SOFC). 

As explained above for a summary or overview, the shaping of all the physically applicable 
properties in a fuel stack system will necessitates agile tools that can integrate multiple physical domains into 
the same model. Equation-based language and causal modeling capabilities enable physics models. This is 
very practical when you want to reuse a component in a variety of settings for a range of objectives. 
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